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ABSTRACT 

The size, mass, luminosity, and space density of Lyman-a emitting (LAE) 
galaxies observed at intermediate to high redshift agree with expectations for 
the properties of galaxies that formed metal-poor halo globular clusters (GCs). 
The low metallicity of these clusters is the result of their formation in low-mass 
galaxies. Metal-poor GCs could enter spiral galaxies along with their dwarf 
galaxy hosts, unlike metal-rich GCs which form in the spirals themselves. Con- 
sidering an initial GC mass larger than the current mass to account for multiple 
stellar populations, and considering the additional clusters that are likely to form 
with massive clusters, we estimate that each GC with a mass today greater than 
2 x 10 5 Mq was likely to have formed among a total stellar mass > 3 x 10 7 M@, 
a molecular mass > 10 9 M , and 10 to 10 9 M & of older stars, depending on the 
relative gas fraction. The star formation rate would have been several M & yr _1 
lasting for ~ 10 7 yrs, and the Lyman-a luminosity would have been > 10 42 erg 
s _1 . Integrating the LAE galaxy luminosity function above this minimum, con- 
sidering the average escape probability for Lya photons (25%), and then dividing 
by the probability that a dwarf galaxy is observed in the LAE phase (0.4%), we 
find agreement between the co-moving space density of LAEs and the average 
space density of metal-poor globular clusters today. The local galaxy WLM, with 
its early starburst and old GC, could be an LAE remnant that did not get into 
a galaxy halo because of its remote location. 



Subject headings: Galaxies: dwarf — Galaxies: star clusters — Galaxies: star 
formation — globular clusters: general 
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Introduction 



Metal-poor globular clusters (GCs) that inhabit the halos of spiral galaxies ( IBrodie fc Strader 



20061 ) could have arrived there in dwarf galaxies t h at got captured by d y namical frictio n 



and dispersed by tidal forces (ISearle fc Zinn 



200C 



same dwarfs ( jGao et al. 



Mackev fc Gilmord 12004; ICarraro et al.l 120071 ; ICasetti-Dinescu. et al. 



1978 



Zinnecker et al 



Many GCs are still in debris streams ( Da Costa &; Armandrofj 1995j ; 



1988L 

iFreemanl 



Palma et al. 



2009; 



1993). 



2002; 



Newberg et al. 



9[~ Mackev et al. 2010) or share orbits with o ther GCs that presumably came from the 



2007: Smith et al. 2009). The low metal abundance of halo GCs 



compared with disk and bulge GCs (IStrader et al.l 120071 ; lAlves-Brito et al.l 1201 ll) c ould then 



be the result of the m ass-met allicity relation for galaxies (IChies-Santos et a 



established early on ( lErb et al.l 120061 ; iMaiolino et al 



metallicity relation is also present for GCs themselves (IWehner et al 



2008 



Mannucci et al 



Mieske et al.ll2010l ). although self-enrichment may be the reason (IBailin fc Harrisll2009 ) . The 



2008 



201l|), which is 
2009T). A mass- 



Forbes et al. 2010 



most heterogeneous GCs could have been dwarf galaxy nuclei (e.g. 
Carretta et aDboioh . 



Bekki fc Norris 2006 



Metal-poor GCs in elliptical galaxies would have formed in the same types of dwarfs as 
metal-p oor GCs in spirals. Som e of these GCs got into the ellipticals after two or more spirals 



merged ( Ashman fc 
formed ( 



Cote et al. 



Zepf 



( IRomanowsky et al 



1998). 



19921 ). while others got in with captured dwarfs after the ellipticals 



Debris streams containing GCs are observed in ellipticals too 
20121 ). T he larger radial distribution of metal-poor GCs compared with 
meta l-rich GCs in ellipticals (jHarrisI 120091 ; iForbes et al.l 1201 ll ; iLiu et al.l 1201 ll ; iFaifer et al. 
201 ll ) presumably results from this difference in their formation sites. 



Metal-rich GCs have a different distribution than metal-poor GCs in spiral galaxies, 
suggesting that they for med in the disks and bulges of these galaxies along with the other 
disk stars (ILarsonl Il988l ). pe rhaps during a clumpy phase (jShapiro et al.l I2010T) when th e 



thick disk and bulge formed (IBeasley et al.ll2002l ; iGriffen et al. 



2010 



Bournaud et all 12009 ) . 



Their higher metallicities result from the generally higher metallicities of their more massive 
hosts. Metal-rich GCs also get into ellipticals during the mergers of spirals, and other metal- 
rich GCs form in ellipticals during the merger itself ( Ashman fc ZepJ 1992 ). Ellipticals then 
have a higher specific frequency of metal-rich GCs than spirals ( IHarrislll99ll ) because of the 
additional GCs that formed during the merger. Ellipticals should also have a higher specific 
frequency of metal-poor GCs if their more central locations in dense galaxy clusters gave 
them better access to a rich population of remnant GCs in cluster dwarfs. 

The mass distri bution functions of metal-poor and metal-rich GCs are about the same 
( IWehner et al.l 120081 ). Both have a peaked shape with the same characteristic mass (~ 
2 x 10 5 M Q ). In the dwarf model for metal-poor GCs, this similarity requires that the com- 
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bined mass distribution of small numbers of clusters from many dwarf galaxies is the same 
as the single mass distribution of many clusters in a large galaxy. Obse rvations confirm this 
simil arity for old massive clusters in local dwarf and spiral galaxies (jKruijssen fc Cooper 
20121 ). The mass function of clusters is ap parently a local property of s tar formation inde- 
pendent of the number of clusters formed ( lElmegreen fc Efremovlll997l ). The change from 
an initial power law or Schechter mass function to a peaked function i s presumably the result 
of selective disruption of low mass clusters shortly after birth (e.g., iParmentier fc Gilmore 



20051: 



Elmegreen 201ol) an d a stellar evaporation rate that is a weak function of cluster mass 



( [McLaughlin fc Fall 



20081 ) 



Here we consider the formation conditions for GCs in modern environments (Sect. 2) 
and assess whether metal-poor GCs are likely to have formed in galaxies that have already 
been observed at high redshift (Sect. 3). The most likely candidates are small active galaxies 
like Lyman a emitters (LAE). Earli er attempts to link GC formation with specific galaxies or 
types of galaxies (e.g.. lLarsonlll988l ) lacked the extensive data on high-redshift star formation 
that are now available in large surveys. A comparison with other models is in Section 4. 



2. Conditions for GC Formation in the Local Universe 



In the solar neighborhood of the Milky Way, star clusters form in the midst of other 
clusters and unbound stars, with a mass that is at most a few percent of the associated 
molecular cloud mass. One of the nearest regions where star formation has produced a total 
stellar mass comparable to that in a globular cluster is Eta Carina, where instead of a single 
cluster there are h undreds of clusters up to ~ 10 3 M & in mass and about three times as much 
mass in free stars (IFeigelson et al.ll201ll ). A bigger star-forming region is W31 at a distance 
of 6 kpc, which emits 6 x 10 6 L & of young star radiation, contains several ultra-compact HII 
regions and subclumps of ~ 10 5 ' 8 L & each, masers and infrared prestellar clumps. The dense 
molecular clumps, however, extend up in mass to only ~ 8000 M & because they are part 
of an dn(M)/dM oc M -1 ' 5 power law of clump masses that distributes the total molecular 
mass of 1.2 x 10 5 M among smaller pieces (Beuther et al. 2011). The largest cluster likely 
to come from such clumps is only several thousand M & . 

Looking for even larger regions in the Milky Way, one of the most extreme is W43, 
located at the near tip of the main stellar bar. Such bar-end regions typically house prodi- 
gious st ar formation, presuma bly because of the massive collection of gas in orbit around 
the bar (IKenney &: Lord! 1 199 ll ). W43 extends for 310 pc and contains a molecular mass of 
7 x 10 6 M in addition to an atomic mass of 3 — 6 x 10 6 M Q (Luong et al. 2011). The 
fraction of molecules in dense 13 CO substructures is high, ~ 10%, and the CO line FWHM 
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is large, 22 km s _1 . This is the type of region that might produce a 10 5 M Q cluster along 
with many smaller clusters (Luong et al. 2011). It took the Milky Way bar to make these 
conditions, which means a total distance for gas collection and compression in the bar flow 
exceeding several kiloparsecs. 



For a power law cluster mass function, dn(M)/dM 
cluster mass to the largest cluster mass is 



n Q 



M P , the ratio of the total 



Mi 



total 



M n 



1 + 



rM M l-f) dM 

M max f™ M-PdM 



for minimum cluster mass M min . This ratio implies that the minimum likely young stellar 
mass to form a single largest cluster of mass M is 



Mi 



total 



r}~ x M + r}~ x {0 ~ 1)M /3_1 



M 2 '■ - Ml 
2-/3 



In 



M 



(2) 



where rj c is the fraction of young stellar mass that forms in bound cluster s as opposed to 
the fie ld, and the second result is for — 2. Observations of Eta Carina by iFeigelson et al. 
( 1201 ll ) suggest r] ~ 0.25. For typical ~ 2 and minimum cluster mass M min ~ 10 M , the 



minimum total young stellar mass to form a 10 6 M Q cluster is 5 x 10 7 M . With an average 
star formation efficiency in molecular clouds of ~ 2%, the associated molecular mass has 
to exceed 2.5 x 10 9 M . We view this as the necessary molecular mass in a whole galaxy, 
since by statistical arguments, subclumps that form individual clusters need not all occur 
in the same shielding envelope for H 2 formation . Similar discu ssio ns leading to a large tota l 
molecular mass for GC formation were made by |L arson! ( 119881 ) and lHarris &: Pudritzl ( 1l994l ). 



For typical gas fractions of 50% or more at high redshift, as determined from CO ob- 
servations or star formation rates convert e d to gas mass using the iKennicuttl ( 119981 ) relation 
( IMannucci et al.l l2009t iDaddi et al.l |2010| ; iTacconi et al.l |2010| ) , this result for gas mass im- 
plies that the stellar mass in a galaxy that forms a single ~ 10 6 M & GC plus all of the other 
clusters and bare stars that accompany it in the starburst i s typically around 10 9 M Q . For 
90% gas fractions at the low metallicities of GCs (Fig. 8 in IMannucci et al.ll2009l ). the un- 
derlying stellar mass might be ~ 10 8 M & . Lower mass clusters would generally be associated 
with lower masses of peripheral stars, and the peripheral mass would be lower still if the GC 
forms in the first galaxy starburst (100% gas fraction). 

An upper mass cutoff for clusters gi ying a Schechter distribution function for mass (i.e., 



a power law with an exponenti al tail - ISchechter 



suggested for spiral galaxies by iGieles et al.l (12 006a 



galaxies by I Waters et al.l ( 120061 ) and I Jordan et al.l ( 120071 ). iLarsenl (120091 ) measured an average 



9761 ) i nstead o f a po wer law has been 



bj) a nd ILarsenl (120091) and for elliptical 
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cutoff (where the exponential tail begins) for several spirals at ~ 2 x 10 5 M . Such a cutoff 
can increase the required cloud mass if the maximum cluster mass in a distribution is larger 
than the cutoff mass. Cutoffs have n ot been explained theoretically except for a possible 
connection with interstellar pressure (jElmegreen &: Elmegreenl l200ll ). The application of 
such a cutoff to the present problem is uncertain bec ause massive young clu sters with M ~ 
10 6 Mp, are observed in dwarf galaxies like NGC 1569 fao fc Filippenkolll996h and NGC 1705 
(IMeurer et al.lll995l ); there is no observation of a statistically significant cutoff in dwarfs. 



The star formation rate for a molecular mass o f 2 x 10 9 is ~ 1 yr 1 for a 



conventional molecular consumption time of 2 Gyrs (jBigiel et al.l 12008c iLeroy et al.ll2008l ). 
This is the characteristic rate for a spiral galaxy today, and large for a dwarf galaxy. If the 
mole cular consumption time was shorter by a factor of ~ 4 in the early universe (IGenzel et al. 



Dekel & Silk 


1986: 


Finkelstein et al. 


2011a 



McLinden et al.ll201lh . 



We are interested here in whether such a star-forming dwarf could be observed at about 
the time when it forms a 10 6 M cluster. For a star formation rate of 4 M Q yr -1 lasting 
~ 100 Myr, which is the orbit time of a 0.5 kpc radius disk with a typical dwarf galaxy 
rotation speed of 30 km s _1 , the SDSS g-band absolute AB magnitude would be —20.0 and 
the young-star luminosity 1 x 10 10 L Q (according to models in Bruzual & Chariot 2003 for 
a metallicity of 0.2 solar; emission lines would make the galaxies brighter). For redshifts 
z = 2 and z = 6, where the distance moduli in the observer frame are 46.0 and 48.9 and the 
restframe g-band shifts to 1.4/im and 3.3/iin, the apparent AB magnitudes would be 26.0 
and 28.9, respectively. These magnitudes are bright enough to be detected in deep surveys. 
Similarly, the absolute g-band AB magnitude of a model cluster 100 Myr old with 10 6 M & 
at birth can be determined from Bruzual & Chariot (2003) to be —12.46. At z = 2 and 6, 
the apparent magnitudes of such a cluster would be 33.5 and 36.4 mag, which are beyond 
the current limits. 

Also related to observability is the diameter of a typical dwarf galaxy, ~ 1 kpc, which 
corresponds to an angular size of 0.12 arcsec at z = 2 and 0.17 arcsec at z — 6. These sizes 
can be resolved by HST and some ground-based instruments with adaptive optics, which 
implies that the formation environments for metal-poor GCs can be observed directly. The 
10 6 M Q clusters themselves would not be resolved, however. 

The co-moving space density of dwarf galaxies that formed GCs during an intense 
starburst should have been about twice the space density of metal-poor GCs today, con- 
sidering one forming GC per dwarf and an evaporative loss of about a factor of 2 dur- 
ing the intervening time (jVesperinil Il998l ). The space density of all GCs today is ~ 8 
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Mpc~ 3 (Portegies Zwart &: McMillan! l2000l ). of which about half to two-thirds are metal- 
poor (IForbes et al.ll2000l ). i.e., the blue GCs in the bimodal color distribution. This makes 
the expected dwarf galaxy co-moving density at the time of metal-poor GC formation equal 
to about 8 Mpc -3 , although only one quarter of these produce GCs more massive than the 
peak in today's GC mass function. We discuss these densities in more detail below. 

The metallicities of metal-poor GCs should be about the same as the gas metallici- 
ties of the gal axies t hey f ormed in, which is low for low-mass galaxies at intermediate to 
high redshift. IZinnl (1198 5) considered that metal poor halo GCs have [Fe/H] < —0.8 and 
Muratov fc Gnedin (2010) considered [Fe/H ] < —1. If we consider the solar abundance to be 
12 + log(0/H) = 8.66 (Usplund et al.lboO^ . then metal-poor GCs have 12 + log(0/H) < 7.7 
or 7.8. This value is consistent with the observed gas metallicities of galaxies with stel- 



lar masses less than 10 9 M and redshifts greater than z ~ 3.5 ( jMaiolino et al.l 12008 



Mannucci et all 120091 ) . 



The local galaxy WLM, at a distance of ~ 1 Mpc ( iMcConnachie et al.l 120051 ). is an 
example of a dwarf with a me tal-poor GC and a very early starburst of about the same age. 
According to iDolphinl ( 120001 ) . the galaxy has an iron abundance that rises from [Fe/H] = 
— 2.18±0.28 to — 1.34 ±0.14 in a prolonged s tarburst that occurred 12 Gyr to 9 Gyr ago, and it 
has a globular clus ter ( lAbles &: Ableslll977l ) with an abundance of [Fe/H]= —1.63 ±0.14 (from 
Hodge et al.l 119990 that also has a very old age, estimated to be ~ 14 Gyr by iHodge et al. 
( 119991 ). Enhanced a eleme nts in the GC suggest a p eriod of rapid star formation in the 
gas from whi ch it formed (IColucci fc Bernstein! 1201 ll ). Most of the WLM halo stars are 
also very old (iMinniti k Ziilstral E997h . The stellar mass of WLM today is ~ 1.6 x 10 7 M© 



( jZhang et al.l |2012[ ); the total dynamical mass within the half-light radius of 1.6 kpc is 
4.3 ± 0.3 x 10 8 Mr,, (|Leaman et al.l 120121). The V - band absolute magnitude of the GC is 
—8.8 mag ( Sandage fc Carlson Il985 ; IHodge et al. 1999 ). which gives it a luminosity 4.8 
times greater than the luminosity at the peak of the Milky Way globular cluster distribution 
function. Sc aling the GC mass ac cordingly makes it 10 6 M , or 10% of the current galaxy 
stellar mass. iLeaman et al.l (120121 ) suggest that the WLM stellar mass was only 8 x 10 5 M Q 
11 Gyr ago, which is comparable to the GC mass at the same time. Evidently, the GC in 
WLM was the dominant part of a major starburst shortly after the galaxy formed. It remains 
inside WLM and not in the halo of th e Milky Way or M31 because of its large distance from 
both, ~ 1 Mpc ( ILeaman et al.l 120121 ) . We propose that most of the other metal-poor GCs 
that formed in young dwarfs in the local group got captured by the Milky Way, M31, or 
M33, or possibly also by the LMC and SMC, before the rest of the dwarf galaxy gas had a 
chance to form more stars like WLM did. Thus the remnant streams may have low stellar 
masses like WLM did when it formed its GC; this is much less than the stellar mass in WLM 
today. 
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In summary, a wide variety of evidence points to the formation of metal-poor globular 
clusters in low mass galaxies at intermediate to high redshift. Simple considerations suggest 
that such galaxies could be observed with modern instruments during their GC-formation 
phase, although the young GCs themselves probably cannot be distinguished. We ask now 
whether the active phases of the host galaxies have been observed already. 



3. Lyman a Emitters as Sites for GC Formation 



Small, metal-poor galaxies at high redshift are a natural choice for the birth place of 
metal-poor globular-clusters. Lyman-a emitting galaxies are typically dwarf star-forming 
galaxies at intermediate to high r edshift. Most LAE galaxies have low stellar masses, typ- 
ically between 10 7 and 10 8 M Jpirzkal et all l2007t iFinkelstein et all bpOTt). although a 
small fraction have higher m ass following a Schechter luminosity function ( jLai et ah 2007 ; 



Finkelstein et al.l 120081 120091 ) . The Lyman a emi tting regions in LAE galaxies are usu 



ally very young, with stellar ages o f a few 10 7 yrs (jPirzkal et al.l 120071 : iGawiser et al.l 12007 



Finkelstein et al 



2007 



2008 



stars too (Nilsson et a 



Pentericci et al. 



2009 



2011 



Bond et al. 



2009 ). Sometimes there is an older underlying population of 

LAEs are usua lly compact (~kpc; 

), with a character- 



Acauaviva et al 



2009 



2010 



2012) 



Finkelstein et al. 



2011 



istic radius that does not evolve much with redshift ( jMalhotra et al.l 120121) . althou gh the 
LAEs with underlying populations could be growing in mass f lNilsson et al.l 120111 ). This 



(Kovac et al.l 


2007 




Gawiser et al. 


2007 




T 


ilvi et al. 


2009 


)• I 


jAEs are also 


poor (~ 0.1 solar; 


Finkelstein et al 


. 2011a; 


Richardson et al. 


2012. in prep) 



LAE galaxy properties can be explained by postulating that most of them are dwarf 
galax ies emitting Lyman-a in the first ~ 3 x 10 7 years after a major starburst (jMalhotra fc Rhoads 
20021 ). The escape of Lyman a radiation would favor a porous interstellar medium in a small 
galaxy, so there may be a selection effect for dwarfs among LAE samples. This selection 
works in our favor because we are looking for dwarfs that had starbursts at intermediate to 
high redshift. We postulate that these dwarfs are the formation sites of metal-poor GCs, 
and the precursors of the remaining stellar streams that accompanied these GCs into galaxy 
halos. Extreme emission line objects at intermediate redshift are in the s ame category, and 



with their enormous starbursts, could have made most of today's dwarfs (Ivan der Wei et al. 
201 lh . 



tion of GCs today is a log-normal distribution centered at ~ 2 x 10 5 M & 
These clusters have lost some of their original mass by evaporation. 



McLaughlin 


2003 




McLaughlin & Fall 
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( 120081 ) fit the current GC mass function with an initial Schechter function and an evaporated 
mass per cluster of A = 1.45 x 10 4 M Q (p h /M Q pc -3 ) 1 / 2 for cluster density p h ~ 10 — 1000 M G 



© 



pc 
A 



inside the half -light radius. For avera ge ph = 246 M pc 3 in the Milky Way, 



2.3 x 10 5 M flMcLaughlin fc Fall 12008 1. The initial GC mass would have been at 
least this much larger than its current mass. 

We estimated above that the number density of original metal-poor clusters had to be 
about - 8 Mpc -3 . This came from a current density of ~ 4 Mpc 3 , multiplied by 2 to 
account for evaporation. For a nearly constant evaporation rate per cluster, this estimate of 
initial density should be done again, now in two parts. First, we consider all clusters with 
a mass today larger than 2 x 10 5 M , which is both the mass at the peak of the GC mass 
function and the total evaporated mass per cluster, on average. An approximation is that all 
clusters originally larger than this are still present today because they have not completely 
evaporated yet. Their density is about half the current density of metal-poor GCs, giving 
~ 2 Mpc -3 , because they represent half of the current GC mass function by number. Thus 
we ask whether galaxies going through the LAE phase and forming clusters more massive 
than 2 x 10 5 M & today had a co- moving space density of around 2 Mpc -3 . 

Second, we consider lower mass clusters. These are on the low-mass, decreasing part 
of the current GC mass function and should be highly evaporated from their initial state. 
The mass function of clusters at birt h is usually observed to be dn(M)/dM oc M -2 (e.g., 



de Grijs &: Anders! 120061 ; lLarsenll2009l ). If we consider as representative the lower part of the 
GC mass function down to 1/4 the mass of the peak, then the initial number of clusters 
below the peak, i.e., between 1/4 and 1 times the current peak, is 3 times the number 
of clusters with masses greater than the current peak, based on this initial cluster mass 
function. Presently, there is an equal number of clusters with a mass below and above the 
peak, so 1/3 of those below the peak survive. The total survival fraction is still the factor 
of 2 estimated in the introduction, but now we have divided it up between 100% survival at 
M > 2 x 10 5 M and 33% survival at initial M < 2 x 10 5 M . (That is, for every GC today 
at M > 2 x 10 5 M there was one GC originally, and for every GC today at M < 2 x 10 5 M 
there were three GCs originally). We first check whether clusters above the current peak 
could have formed in LAE galaxies. 

To make a cluster of mass 2 x 10 5 M today, the initial stellar mass had to be larger, 
not only because of evaporation over a Hubble time, as discussed above, but also to account 
for the high abundance of elements from first-generation stellar ejecta that are present in 
second- generation stars. Globular c lusters have multiple stellar populations ( see rev i ews i n 
Carretta et aDboiObh . Models by b'Ercole et aD (booa . koii ). iBekkil fcpllh 



Rcnzini 



2008 



Schaerer &: CharbonnellEoil ) and others suggest that the initial cluster mass was ~ 10 6 M 
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or more, ten times larger than the present day mass. These models are uncertain, but to be 
consistent with them, we assume that the initial characteristic cluster had 2 x 10 6 M . It 
would most likely have been accompanied by other clusters, making the total young stellar 



mass ~ 3 x 10 7 M & if r\ — 1 in equation 



Most of the star formation that made these 



clusters, whether in the first or the second generations, wo uld presumably occu r in ~ 10 7 
years, before supernovae could remove the star- forming gas (ID'Ercole et al.ll2012l ). However, 
the total span of cluster formation could have taken ~ 100 Myr, considering the need for 
enriched ejecta from first-generation stars to re-accumulate along with fresh material to form 
the second generation. Non-clustered stars would most likely have been present with the 
first generation too (77 < 1), but they would presumably be relatively dim by the time the 
second generation stars formed in the final cluster burst. Considering the many uncertainties 
that are still present in the detailed theory of GC formation, we use here a representative 
luminosity for the starburst phase that is equivalent to 3 x 10 7 M & of stars forming in ~ 10 
Myr. This gives a star- formation rate during the main starburst equal to 3 M yr _1 . 

Another consideration is the absorption by dust of Lyman continuum radiation before 
it ionizes the nearby gas and produces Lyman a radiation. We assume a small amount, 
Air, y 1 magnitude of ex tinction for Lyman continuum radiation (e.g., iFinkelstein et al. 



201ld ; [Walter et al.l 120121 ) . That reduces the effective star formation rate for Lyman a 
production to ~ 1 M Q yr _1 . This implies a Lyman-a line luminosity of 1 x 10 42 erg s -1 if the 
escape fraction of the Lyman-a photons is 100%. Recent studies of LAE galaxies show that 



the average escape fraction of Lyman-o; photons is between 17% and 30% (IBlanc et al.ll2011 



Zheng et al.ll2012t iRichardson et al.ll2012. in prepl ). Assuming an average escape fraction of 
25%, a SFR of 1 M /year translates into Lyman-a line luminosity of 2.5 x 10 41 ergs -1 . 



Now we integrate the Lyman-a luminosity function (IZheng et al.ll2012bl ) with log(L*) = 
42.86 ± 0.06 and 0* = —3.55 ± 0.09 down to this minimum luminosity of 2.5 x 10 41 ergs _1 . 
We assume a slope of —2 at low mass. This gives the number density of Lyman-o: galaxies 
brighter than the cutoff for forming a characteristic GC equal to 7 x 10~ 3 Mpc~ 3 . 

In deriving this density, we have assumed that the Lyman-a phase lasts for 10 7 years, 
whereas we observe that the Lyman-a luminosity function is essentially unchanged from 
z = 2 to 7, a time interval of 2.56 Gyrs. This duration for the LAE morphology corresponds 
to ~ 256 generations of individual Lyman-a events, or a probability of catching a low-mass 
galaxy in the LAE phase equal to 0.39%. Dividing the observed number density of LAE 
galaxies that can make today's M > 2 x 10 5 M & GCs by the probability of observation 
produces a total number density of such galaxies that go through the LAE phase equal to ~ 
1.9 Mpc -3 . If each LAE phase forms a single GC with a present-day mass of M > 2 x 10 5 M Q 
(along with smaller clusters, as discussed above), then this would be the comoving number 
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density of metal-poor GCs that end up more massive than the median mass of today's GCs. 
This is in fact the initial co-moving density of low-metal, high-mass GCs that we expected 
from the above discussion of current GC density and evaporation rate. Thus dwarf LAE 
galaxies at redshifts from 2 to 7 could have formed most of the metal poor GCs observed 
today in the high-mass half of the GC mass distribution function. 

Formation of the low-mass GCs requires a higher space density of lower-luminosity LAEs 
because there were three time as many initially low-mass GCs as high mass GCs. This result 
is easily obtained, however, because the LAE luminosity function, presumed equal to the 
LAE mass function, has the same slope as the initial GC mass function, namely, —2 for 
equal intervals of mass. Thus there were 3 times as many LAEs between 1/4 and 1 times 
the fiducial luminosity estimated above as there were above the fiducial luminosity, just as 
there were 3 times as many initial GCs below the current peak as above. We conclude that 
all of the GCs could have been made in LAE galaxies, based these assumptions. 

The above estimates can be summarized by two equations: 

n(M > M %min ) = / $(L)dL = $T(a + 1, L mm /L*) (3) 

^LyA JL 

ram 

and 

T _ M^ min in -o.4A LC 10 42 erg s' 1 

i-'min — , J esc,LyA Ly J A 7T i ; y±) 

tLyA ■ M Q yr 1 

where M* im j n p» 13M c i from equation fl2]) is the minimum stellar mass formed during the 
starburst lasting for the time t^ y A- We assume a luminosity function described by a Schechter 
function. Our calculation above used fiducial values of T = 2.56 Gyrs, t Ly A = 10 7 yr, 
fescLyA ~ 0.25, A LC ~ 1, log($* x cMpc 3 ) « -3.55, log^/ergs" 1 ) « 42.86, a = -2, and 
M c i « 2 x 1O 6 M . 

These equations show that the space density of GCs derived from the LAE luminosity 
function is insensitive to the lifetime tj^A of the Lyman-a phase. If we assume that the life 
span of the Lyman-a emitting phase goes up to t^yA = 3 x 10 7 years, then the minimum SFR 
and the Lyman-a luminosity required would be lower by a factor of three. The number of 
Lyman-a emitters brighter than this cutoff would increase by roughly a factor of 3. However, 
the total number of GCs produced over the redshift interval z = 2 to 7 would remain the 
same, since our 2.56 Gyr interval is now only ~ 85 generations of galaxies. 



Formal uncertainties in the luminosity function parameters taken from IZheng et al. 



( 12012bl ) change the final numbers by only about 6%, since the errors in 0* and L* are 
anticorrelated and roughly cancel out. There is however a larger spread in the luminosity 
function parameter found for LAEs by different authors and surveys. Taking two extremes: 
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(1) iHu et al.l (120101 ) found 0* = —3.96 and L* = 43.0 at z=5.7, which would yield about half 
of the GC number density calculated above; (2) ICiardullo et al.l (120121 ) found 0* = —2.96 and 
L* = 42.6 at z=3.1 yielding 2.3 times the number density derived above. A change of slope 
in the faint end of the luminosity function from a = —2.0 to a = —1.8 reduces the number 
of GCs to 60% of the ab ove value. The fai nt end of the luminosity func tion is quite steep, 
however, in both LAEs rtRauch et all l2008h and Lyman break galaxies (lYan fc Windhorstl 
20041 ); thus, —2 may be a more reasonable estimate for the slope. 



The Lyman-a escape fraction is also uncertain, perhaps by a factor of 1.6 to 2, as is the 
extinction for Lyman Continuum radiation, A-^c, the initial representative cluster mass, M c \, 
and the fraction of stars formed in clusters, rj. A lower escape fraction, higher A-^c, lower 
M c i, or higher 77 will all lower the observed Lyman-a luminosity from a GC-forming galaxy, 
and this increases the space density of such galaxies according to the luminosity function. 
Th e resulting d e nsity of GCs increases in proportion. Lower M c \, as recently proposed 
by lLarsen et al.l (120121 ). can be accompanied by a proportionally lower cluster formation 
fraction, rj, and not change the GC space density much. The formation of more than one 
GC per LAE galaxy would not increase the space density of GCs, because we are accounting 
already for all of the luminosity associated with GC formation. If 2 GCs are made in each 
LAE, then the LAEs would have twice the fiducial luminosity (or each would form a single 
GC twice as often), and about half the estimated number of LAEs would be forming these 
GCs. The total number of final GCs remains about the same. 

Other small galaxies in the early universe could be formation sites for GCs too, but 
if these galaxies are not Lyman-a emitters, then they are probably not in the starburst 
phase at the time of observation, and whatever GCs they formed came about earlier (or 
may come about later). In that case, they are already counted in the above discussion 
when we considered the ~ 0.39% probability for the LAE phase. In addition, some dwarf 
starburst galaxies th at form metal-poor GCs may have low Lyman a escape fractions (e.g. 
Heckman et al.l 1201 ll ). placing them below our limit of integration over the luminosity func- 
tion. These galaxies balance those with high escape fractions, since we use only the average 
value of 25% here as an estimate for this fraction. More massive galaxies would have formed 
GCs in starbursts as well, but then the GC metallicities would be larger at the redshifts 
when typical galaxies became massive, and the result would be metal-rich GCs. This is the 
likely formation mechanism for met al-rich GCs, which a re still correlated with the disk and 
bulge populations of spiral galaxies (jShapiro et al.ll2010l ). 



4. Comparison with other Models 



P revious models of GC formation in a cosmological context (e.g 



2002 



Beaslev et al 



Griffen et al. 



2010 



2002 



Kravstov fe Gnedin 2005: Moore et al. 2006: Bekki et al. 



Bromm & ( 


Ularke 


Bekki et al. 


2008; 



Muratov fe Gnedinll2010[ ) contain many aspects of the present proposal 



witho ut identifying likely hosts that could be visible at intermediate to high redshift. iBeasley et al 
( 120021 ) discussed metal-poor GC formation in gassy disks at early times and suggested that 
the GC formation sites might be damped Lya systems, which are gas clou ds seen in ab- 



sorpt ion with high column densities and very low emission from young stars (IRafelski et al. 
201ll ). We suggest that GC formation sites are brighter than this. 



A variety of models for t riggering GC formation h ave been considered, such as colli- 
sions between dwarf galaxi es (IMuratov fe Gnedin! |2010| ) or pressures from a nearby galaxy 
( IGray fe Scannapiecoll201ll ). Here we discussed GC formation in small galaxies regardless of 
their proximity to large galaxies and regardless of what triggers the star formation. If LAE 
galaxies are the formation sites, then any correlation between LAEs and large galaxies in the 
early universe could reveal the extent to which metal-poor GCs formed in halo sub-clumps 
versus the field. The observation of isolated LAEs would imply that their metal-poor GCs 
entered today's galaxy halos later, presumably making the observed remnant stellar streams 
by tidal disruption. 

Observations show overdensities of LAE galaxies only on the scale of ga l axy clusters. 



namely, ~ 10 and several Mpc flKurk et al 



2007: Gawiser et al. 2007; Tilvi et al. 2009:lQverzier et al 



2004 



2006 



Venemans et al. 2005; Kovac et al. 



20081 ) . This cluster overden- 



sity allows centrally-located elliptical galaxies to accrete a higher proportion of metal-poor 
GCs from former LAEs than spiral galaxies can. Such enhanced accretion is consistent with 
a greater specific frequency of metal-poor GCs in ellipticals compared to spirals. 

For dwarf galaxy hosts in general, metal-poor GCs should form over a range of times 
because their abundance reflects the metal-poor env ironment of the dwarf galaxy, r ather than 
the a yerage metallicity of th e universe. NGC 1569 (jKobumicky fe Skillmanlll9971 ) and NGC 
1705 ( ILee fe Skillmanl 12004 ) are examples of moderately metal-poor environments that make 
massive clusters today. They did not fall into the halos of more massive spiral galaxies, and 
so their stellar masses and heavy element abundances continued to increase. Other models 
have been concerned with truncating the f ormation of GCs a fter a short formation epoch 
in order to keep the metallicity low (e.g., iBekki et al.l 120081 ). A key observation will be 



the age spread of metal-poor GCs. Extended format i on tim es for metal-poor GCs were 
recently observed for SO galaxies by IChies-Santos et al.l (1201 lh . There was als o a suggestion 



of ext ended metal-poor GC formation in one of the stellar population fits by iMendel et al. 
j2007h for Milky Way GCs. 
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5. Conclusions 

Metal-poor GCs could have entered the halos of spiral galaxies in the form of starburst 
remnants in old dwarf galaxies. Their low metallicities are not the exclusive result of an 
earlier birth time compared to metal-rich disk and bulge GCs, but rather the result of a 
lower mass host, considering the mass-metallicity relation in galaxies. Based on observations 
of Milky Way star-formation and gas fractions for high redshift galaxies, we propose that a 
typical metal-poor GC with an original mass of 2 x 10 6 M formed among a total cluster 
mass of ~ 3 x 10 7 M Q along with several xlO 9 M of dense gas and ~ 10 7 to ~ 10 9 M of 
older stars. During a ~ 100 Myr burst phase when the star formation rate was several M Q 
yr _1 , such galaxies would have AB magnitudes of ~ 26 and ~ 29 at wavelengths of 1.4//m 
and 3.3/zm for z ~ 2 and z — 6, respectively, and the GCs would have magnitudes of ~ 34 to 
~ 36 in this redshift range. The host dwarf galaxies would be ~ 1 kpc in size and subtend 
an angle of ~ 0.15". 

Lyman-a emitting galaxies are plausible formation sites for metal-poor GCs. They have 
the right size, luminosity, star formation rate, and metallicity. To determine whether LAEs 
also have the right space density to form all of today's metal-poor GCs, we integrated the 
LAE luminosity function above ~ 2.5 x 10 41 ergs s _1 , which is the expected Lyman a escape 
luminosity of a star-forming region that produces a typical GC today. This limit assumes an 
average 25% Lyman-a escape fraction and one magnitude of extinction for Lyman continuum 
radiation. The result equals the space density of today's metal-poor GCs when corrected for 
the fraction of LAE galaxies that are in the active phase (~ 0.4%). 

JER and SM acknowledge support from NSF AST-0808165. Helpful comments by the 
referee are gratefully acknowledged. 

REFERENCES 

Abies, H.D. & Abies, P.G. 1997, ApJS, 34, 245 

Acquaviva, V., Vargas, C, Gawiser, E., & Guaita, L. 2012, ApJL, 751, 26 

Alves-Brito, A., Hau, G.K.T., Forbes, D.A., Spitler, L.R., Strader, J., Brodie, J. P., Rhode, 
K.L. 2011, MNRAS, 417, 1823 

Ashman, K.M., Zepf, S.E. 1992, ApJ, 384, 50 

Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C, & Kiselman, D. 2004, A&A, 
417, 751 



- 14 - 



Bailin, J., & Harris, W.E. 2009, ApJ, 695, 1082 

Beasley M. A., Baugh C. M., Forbes D.A., Sharpies R. M., Frenk C. S., 2002, MNRAS, 333, 
383 

Bekki K., & Norris J. E., 2006, ApJL, 637, 109 

Bekki K., Yahagi H., Nagashima M., & Forbes D.A., 2008, MNRAS, 387, 1131 
Bekki, K. 2011, MNRAS, 412, 2241 

Bigiel, F., Leroy, A., Walter, F., Brinks, E., de Blok, W. J. G., Madore, B., & Thornley, M. 
D. 2008, AJ, 136, 2846 

Blanc, G. A., et al. 2011, ApJ, 736, 31 

Bond, N. A., Feldmeier, J. J., Matkovic, A., Gronwall, C, Ciardullo, R., & Gawiser, E. 2010, 
ApJ, 716, L200 

Bond, N. A., Gawiser, E., Gronwall, C, Ciardullo, R., Altmann, M., & Schawinski, K. 2009, 
ApJ, 705, 639 

Bournaud, F., Elmegreen, B.G., & Martig, M. 2009, ApJL, 707, 1 
Brodic, J. P., & Strader, J. 2006, ARA&A, 44, 193 
Bromm, V., & Clarke, C. J. 2002, ApJ, 566, LI 
Bruzual, C, & Chariot, S. 2003, MNRAS, 344, 1000 
Carraro, G., Zinn, R., & Moni Bidin, C. 2007, A&A, 466, 181 

Carretta, E., Bragaglia, A., Gratton, R.G., Lucatello, S., Bellazzini, M., Catanzaro, G., 
Leone, F., Momany, Y., Piotto, G., & D'Orazi, V. 2010, ApJ, 714, L7 

Carretta E., Bragaglia A., Gratton, R.G., Recio-Blanco, A., Lucatello S., DOrazi V., & 
Cassisi S., 2010, A&A, 516, 55 

Casetti-Dinescu, D.I., Girard, T.M., Majewski, S.R., Vivas, A. K., Wilhelm, R., Carlin, J.L., 
Beers, T.C., van Altena, W.F. 2009, ApJL, 701, 29 

Chies-Santos, A.L., Larsen, S.S., Kuntschner, H., Anders, P., Wehner, E.M., Strader, J., 
Brodie, J.P, & Santos, Jr., J.F.C. 2011, A&A, 525, A20 

Ciardullo, R., Gronwall, C, Wolf, C, et al. 2012, ApJ, 744, 110 



- 15 - 

Cohen, J. G., Blakeslee, J. P., & Ryzhov, A. 1998, ApJ, 496, 808 

Colucci, J.E., & Bernstein, R.A. 2011, EAS Publications Series 48, (Cambridge: Cambridge 
Univ. Press), 275 

Cote, P., Marzke, R.O., & West, M.J. 1998, ApJ, 501, 554 

Da Costa, G. S., Armandroff, T. E. 1995, A J, 109, 2533 

Daddi, E., Bournaud, F., Walter, F., et al. 2010, ApJ, 713, 686 

de Grijs, R. & Anders, P. 2006, MNRAS, 366, 295 

DErcole A., Vesperini, E., DAntona, F., McMillan, S.L.W., & Recchi, S., 2008, MNRAS, 
391, 825 

DErcole A., DAntona F., Carini, R., Vesperini, E., & Ventura, P. 2012, MNRAS, 423, 1521 

Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

Dolphin, A.E. 2000, ApJ, 531, 804 

Elmegreen, B.G. & Efremov, Y.N. 1997, ApJ, 480 235 

Elmegreen, B.G., & Elmegreen, D.M. 2001, AJ, 121, 1507 

Elmegreen, B.G. 2010, ApJ, 712, L184 

Erb D. K., Shapley A. E., PettiniM., Steidel, C.C., Reddy, N.A., Adelberger, K.L., 2006, 
ApJ, 644, 813 

Faifer, F.R., Forte, J.C., Norris, M.A., Bridges, T., Forbes, D.A., Zepf, S.E., Beasley, M., 
Gebhardt, K., Hanes, D.A., Sharpies, R.M. 2011, MNRAS, 416, 155 

Feigelson, E.D. et al. 2011, ApJS, 194, 9 

Finkelstein, S. L., Rhoads, J. E., Malhotra, S., & Grogin, N. 2009, ApJ, 691, 465 
Finkelstein, S. L., Rhoads, J. E., Malhotra, S., Grogin, N., & Wang, J. 2008, ApJ, 678, 655 
Finkelstein, S. L., Rhoads, J. E., Malhotra, S., Pirzkal, N., & Wang, J. 2007, ApJ, 660, 1023 
Finkelstein, S. L., et al. 2011a, ApJ, 729, 140 

Finkelstein, S.L., Cohen, S.H., Windhorst, R.A. et al. 2011b, ApJ, 735, 5 



- 16 - 

Finkelstein, S.L., Cohen, S.H., Moustakas, J., Malhotra, S., Rhoads, J.E., & Papovich, C. 
2011, ApJ, 733, 117 

Forbes, D.A., Masters, K.L., Minniti, D., & Barmby, P. 2000, A&A, 358, 471 

Forbes, D.A., Spitler, L.R., Harris, W. E., Bailin, J., Strader, J., Brodie, J. P., Larsen, S.S. 
2010, MNRAS, 403, 429 

Forbes D., Spitler L., Strader J., Romanowsky A., Brodie J., Foster C, 2011, MNRAS, 413, 
2943 

Freeman, K.C. 1993, in The Globular Cluster- Galaxy Connection, ASP Conference Series 
48, eds. G.H. Smith and J.P Brodie, 608 

Gao, S., Jiang, B.-W., & Zhao, Y.-H. 2007, ChJAA, 7, 111 

Gawiser, E., et al. 2007, ApJ, 671, 278 

Genzel, R., Tacconi, L. J., Gracia-Carpio, J., et al. 2010, MNRAS, 407, 2091 
Gieles, M., Larsen, S.S., Bastian, N., & Stein, I.T. 2006a, A&A, 450, 129 
Gieles, M., Larsen, S.S., & Sheepmaker, R.A. 2006b, A&A, 446, L9 
Gray, W.J., & Scannapieco, E. 2011, ApJ, 742, 100 

Griffen, B. F., Drinkwater, M. J., Thomas, PA., Helly, J.C., & Pimbblet, K. A. 2010, 
MNRAS, 405, 375 

Harris, W.E. 1991, ARAA, 29, 543 

Harris W., 2009, ApJ, 699, 254 

Harris, W.E., & Pudritz, R.E. 1994, ApJ, 429, 177 

Heckman, T.M., Borthakur, S., Overzier, R., Kauffmann, G., Basu-Zych, A., Leitherer, C, 
Sembach, K., Martin, D.C., Rich, R.M., Schiminovich, D., & Seibert, M. 2011, ApJ, 
730, 5 

Ho, L.C., & Filippenko, A.V. 1996, ApJ, 472, 600 

Hodge, P.W., Dolphin, A.E., Smith, T.R., & Mateo, M. 1999, ApJ, 521, 577 

Hu, E. M., Cowie, L. L., Barger, A. J., Capak, P., Kakazu, Y., Trouille, L. 2010, ApJ, 725, 
394 



-17- 



Jordan, A., McLaughlin, D.E., Cote, P., Ferrarese, L., Peng, E.W., Mei, S., Villegas, D., 
Merritt, D., Tonry, J.L., & West, M.J. 2007, ApJS, 171, 101 

Kenney, J.D.P, & Lord, S.D. 1991, ApJ, 381, 118 

Kennicutt, R.C. 1998, ApJ, 498, 541 

Kobulnicky, H.A., & Skillman, E.D. 1997, ApJ, 489, 636 

Kovac, K., Somerville, R. S., Rhoads, J. E., Malhotra, S., & Wang, J. 2007, ApJ, 668, 15 
Kravstov A. V. & Gnedin O. Y., 2005, ApJ, 623, 650 
Kruijssen, J.M.D., & Cooper, A.P 2012, MNRAS, 420, 340 

Kurk, J. D., Pentericci, L., Rdttgering, H. J. A. & Miley, G. K. 2004, A&A, 428, 793 
Lada, C.J., & Lada, E.A. 2003, ARAA, 41, 57 

Lai, K., Huang, J.-S., Fazio, G., Cowie, L. L., Hu, E. M., & Kakazu, Y. 2007, ApJ, 655, 704 
Larsen, S.S. 2009, A&A, 494, 539 

Larsen, S.S., Strader, J., & Brodie, J.P 2012, A&A, in press, arXiv 1207.5792 

Larson, R.B. 1988, in The Harlow-Shapley Symposium on Globular Cluster Systems in 
Galaxies, eds. J.E. Grindlay & A.G.D. Philip, IAUS 126, (Dordrecht, Kluwer), 311 

Leaman, R., Venn, K.A., Brooks, A.M., Battaglia, G., Cole, A. A., Ibata, R.A., Irwin, M.J., 
McConnachie, A.W., Mendel, J.T., & Tolstoy, E. 2012, ApJ, 750, 33 

Lee, H., & Skillman, E. D. 2004, ApJ, 614, 698 

Leroy, A.K., Walter, F., Brinks, E., Bigiel, F., de Blok, W. J. G., Madore, B., & Thornley, 
M. D. 2008, AJ, 136, 2782 

Liu C, Peng E., Jordan A., Ferrarese L., Blakeslee J., Cote, P., Mei S., 2011, ApJ, 728, 116 

Mackey A. D., & Gilmore G. F., 2004, MNRAS, 355, 504 

Mackey, A. D., Huxor, A. P., Ferguson, A.M.N., Irwin, M.J., Tanvir, N. R., McConnachie, 
A. W., Ibata, R. A., Chapman, S. C, Lewis, G. F., 2010, ApJ, 717, Lll 

Maiolino R. et al, 2008, A&A, 488, 463 

Malhotra, S., & Rhoads, J. E. 2002, ApJ, 565, L71 



Malhotra, S., Rhoads, J. E., Finkelstein, S. L., Hathi, N., Nilsson, K., McLinden, E., k 
Pirzkal, N. 2012, ApJ, 750, L36 

Mannucci F. et al, 2009, MNRAS, 398, 1915 

McConnachie, A. W., Irwin, M. J., Ferguson, A. M. N., Ibata, R. A., Lewis, G. F., Tanvir, 
N. 2005, MNRAS, 356, 979 

McLaughlin, D.E. 2003, in Extragalactic Globular Cluster Systems, ESO Astrophysics Sym- 
posia, ed. M. Kissler-Patig. Berlin: Springer- Verlag, p. 329 

McLaughlin, D. E., & Fall, S. M. 2008, ApJ, 679, 1272 

McLinden, E.M., Finkelstein, S.L., Rhoads, J.E. et al. 2011, ApJ, 730, 136 

Mendel, J.T., Proctor, R.N., k Forbes, D.A. 2007, MNRAS, 379, 1618 

Meurer, G.R., Freeman, K. C, Dopita, M.A., k Cacciari, C. 1992, AJ, 103, 60 

Meurer, G. R., Heckman, T. M., Leitherer, C, Kinney, A., Robert, C. k Garnett, D. R. 
1995, AJ, 110, 2665 

Mieske, S. et al. 2010, ApJ, 710, 1672 

Minniti, D., k Zijlstra, A.A. 1997, AJ, 114, 147 

Moore B., Diemand J., Madau P., Zemp M., k Stadel J. 2006, MNRAS, 368, 563 
Muratov, A.L., k Gnedin, O.Y. 2010, ApJ, 718, 1266 
Nakasato, N., Mori, M., k Nomoto, K. 2000, ApJ, 535, 776 
Newberg, H.J., Yanny, B., k Willett, B.A. 2009, ApJ, 700, L61 

Nilsson, K.K., Ostlin, G., M o Her, P., Moller-Nilsson, O., Tapken, C, Freudling, W. k 
Fynbo, J.PU 2011, A&A 529, A9 

Overzier, R.A., Miley, G.C., Bouwens, R.J. et al. 2006, ApJ, 637, 58 

Overzier, R.A., Bouwens, R.J., Cross, N. J. G., et al. 2008, ApJ, 673, 143 

Palma C, Majewski S. R., k Johnston K.V., 2002, ApJ, 564, 736 

Parmentier, G. k Gilmore, G. 2005, MNRAS, 363, 326 



- 19 - 



Pentericci, L., Grazian, A., Fontana, A., Castellano, M., Giallongo, E., Salimbeni, S., & 
Santini, P. 2009, A&A, 494, 553 

Pirzkal, N., Malhotra, S., Rhoads, J. E., & Xu, C. 2007, ApJ, 667, 49 

Portegies Zwart, S. F., & McMillan, S. L. W. 2000, ApJ, 528, L17 

Rafelski, M., Wolfe, A.M., & Chen, H.-W. 2011, ApJ, 736, 48 

Rauch, M. et al. 2008, ApJ, 681, 856 

Renzini A., 2008, MNRAS, 391, 354 

Richardson, M. L. A., McLinden, E. M., Malhotra, S., Rhoads, J. E., & Hibon, P. 2012, in 
prep, ApJ 

Romanowsky, Aaron J.; Strader, Jay; Brodie, Jean P.; Mihos, J. Christopher; Spitler, Lee 
R.; Forbes, Duncan A.; Foster, Caroline; Arnold, Jacob A. 2012, ApJ, Journal, 748, 
29 

Sandage, A., & Carlson, G. 1985, AJ, 90, 1464 
Schaerer, D., & Charbonnel, C. 2011, MNRAS, 413, 2297 
Schechter, P. 1976, ApJ, 203, 297 
Searle L., & Zinn R. 1978, ApJ, 225, 357 

Shapiro, K.L., Genzel, R., & Fdrster Schreiber, N.M. 2010, MNRAS, 403, L36 

Smith, M. C, Evans, N. W., Belokurov, V., Hewett, P.C., Bramich, D. M., Gilmore, G., 
Irwin, M. J., Vidrih, S., & Zucker, D. B. 2009, MNRAS, 399, 1223 

Strader, J., Beasley, M.A., & Brodie, J.P 2007, AJ, 133, 2015 

Tacconi, L. J., Genzel, R., Neri, R., et al. 2010, Nature, 463, 781 

Tilvi, V., Malhotra, S., Rhoads, J. E., Scannapieco, E., Thacker, R. J., Iliev, I. T., & 
Mellema, G. 2009, ApJ, 704, 724 

van der Wei, A. et al. 2011, ApJ, 742, 111 

Venemans, B. P., Rottgering, H.J. A., Miley, G.K., Kurk, J.D., De Breuck, C, Overzier, 
R. A., van Breugel, W. J. M., Carilli, C. L., Ford, H., Heckman, T., Pentericci, L., 
McCarthy, P. 2005, A&A, 431, 793 



-20 - 

Vesperini, E. 1998, MNRAS, 299, 1019 

Walter, F., Decarli, R., Carilli, C, Riechers, D., Bertoldi, F., Weiss, A., Cox, P., Neri, R., 
Maiolino, R., & Ouchi, M. Egami, E., & Nakanishi, K. 2012, ApJ, 752, 9 

Waters, C.Z., Zepf, S.E., Lauer, T.R., Baltz, E.A., & Silk, J. 2006, ApJ, 650, 885 

Wehner, E.M.H., Harris, W.E., Whitmore, B.C., Rothberg, B., & Woodley, K.A. 2008, ApJ, 
681, 1233 

Yan, H. & Windhorst, R.A. 2004, ApJ, 600, LI 

Zhang, H.-X., Hunter, D. A., Elmegreen, B. C, Gao, Y., & Schruba, A. 2012, AJ, 143, 47 
Zheng, Z.-Y., et al. 2012b, ArXiv 1111.3386 
-. 2012, ApJ, 746, 28 
Zinn, R. 1985, ApJ, 293, 424 

Zinnecker, H., Keable, C. J., Dunlop, J.S., Cannon, R.D., Griffiths, W. K. 1988, IAUS, 126, 
603 



This preprint was prepared with the A AS IATgX macros v5.2. 



